function [co] = convert_rhombohedral_to_hexagonal(co)
%CONVERT_RHOMBOHEDRAL_TO_HEXAGONAL Summary of this function goes here
%   Detailed explanation goes here


A_r = co.LatticeVectors(1,:);
B_r = co.LatticeVectors(2,:);
C_r = co.LatticeVectors(3,:);

A_h = A_r - B_r;
B_h = B_r - C_r;
C_h = A_r + B_r + C_r;

% A_r = [0 norm(A_h)*sqrt(3) norm(C_h)]/2;
% B_r = [0 norm(A_h)*sqrt(3) norm(C_h)]/2;
% C_r = [0 norm(A_h)*sqrt(3) norm(C_h)]/2;

co.LatticeVectors = [A_h; B_h; C_h];

end



